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Abstract 

We compute analytically the complete electroweak two-loop corrections to the B° — B° 
mixing. These corrections fix the normalization of the electroweak coupling employed in 
the extraction of \V t d\ and reduce the theoretical uncertainty due to higher order elec- 
troweak effects from several percent to a few parts in a thousand. If the LO result is 
expressed in terms of or of the MS coupling g(M z ), the two-loop corrections are 
0(1%), the exact value depending on the mass of the Higgs boson. We discuss in detail 
the renormalization procedure and the scheme and scale dependence, and provide practi- 
cal formulas for the numerical implementation of our results. We also consider the heavy 
top mass expansion and show that in the case at hand it converges very slowly. 



1 Supported by the Bundcsministcrium fur Bildung und Forschung under contract 06 TM 874 and by 
the DFG project Li 519/2-2. 



1 Introduction 



The B° — B° system offers rich possibilities for studies of CP violation and the quark 
mixing structure of the Standard Model (see Jj]] for a recent review). The physics of this 
system is well described by Ti.^ =2 , the effective low-energy Hamiltonian for the B° <-> B° 
transition, and the most important observable directly linked to Ti,^ =2 is AM B °, the 
mass difference between the heavy and the light mass eigenstates in the B° — B° system. 
Theoretically, this quantity is given by 

AM B o = —\(B°\nfff 2 \B°)\, (1) 
rriB 

and the experimentally measured value is || 

AM B o = (0.46 ± 0.02) x 10 12 s -1 . (2) 

If all other ingredients in the evaluation of the r.h.s. of Eq. (]!]) are sufficiently well known, 
one can extract from this measurement the absolute value of the Cabibbo-Kobayashi- 
Maskawa (CKM) parameter Vtd, which plays an important role in the standard analysis 
of the unitary triangle [3J. This example illustrates the phenomenological relevance of 
a precise determination of Ti,^ =2 and its matrix element, especially if one takes into 
account the experimental accuracy already achieved in the measurement of AM B o ■ 

The two main steps for such a determination are (i) the calculation of Ti^ =2 (or more 
precisely: the corresponding Wilson coefficient) in perturbation theory - in particular, this 
calculation yields the complete dependence of H^ =2 on all the heavy degrees of freedom, 
and (ii) the evaluation of the remaining low-energy matrix element on the lattice or 
by some other non-perturbative method. Presently, this second step causes the largest 
theoretical uncertainty in the determination of (Ti,^- =2 ), about ±20% (dominated by 
systematics of the lattice calculation, i. e. extrapolation to the continuum ||,(4]]). There is, 
however, well-founded hope that with increasing computer power or by new developments 
in lattice theory, this uncertainty may be significantly reduced in the near future. 

The perturbative calculation, on the other hand, has been by now performed at the 
level of next-to-leading order QCD ||. This analysis includes the matching at 0(a s ) 
and the 0(a 2 ) anomalous dimension of the relevant four-quark operator; the achieved 
accuracy is better than ±1% ||. 

In this paper, we focus on a different part of the perturbative analysis which has not 
yet been considered in the literature: the two-loop electroweak correction to the effective 
Hamiltonian Ti,^ =2 . In our opinion, at least three reasons suggest such a calculation: 

• Reduction of scheme dependence. The leading order (LO) result (Fig. la shows one of 
the relevant diagrams) is proportional to g 4 , where g is the weak coupling. However, 
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Figure 1: (a) shows one of the box diagrams contributing to the B° — > B° transition at 
leading order, whereas (b) depicts an example of an electroweak two-loop correction to this 
process which may be seizable in the large M t limit. 

at LO no renormalization prescription for g is needed, so one can use the numerical 
value of g in any renormalization scheme. For example, the difference between g A 
at the scale M z calculated from sin 2 ^ = 0.23155 and a(M z ) = (128.9)- 1 f§, 
or from the relation g 2 = 8M^G^/^/2 amounts to about 2.5%. As noted in M, 
such ambiguity reflects the uncertainty of the LO result due to the uncalculated 
electroweak two-loop correction, and it already exceeds the existing perturbative 
QCD uncertainty mentioned above. In a similar way, if we express the LO result in 
terms of the MS coupling #(//), and change the scale \i between M w /2 and 2M W , 
we obtain a 5% variation, a normalization ambiguity which is almost completely 
removed by the consideration of two-loop electroweak effects. 

• Possible large corrections due to a heavy top. The coupling of the Higgs particle to 
the top quark is proportional to gM t /M w . Since gM t /M w « g s at the scale M z , 
from diagrams like the one in Fig. lb one can expect large contributions of roughly 
the same order of magnitude as the QCD corrections. 

• Performing a complete two-loop electroweak calculation. Despite recent efforts, there 
exist very few nearly complete electroweak calculations at the two-loop level flOfl , 
and many available results rely on heavy mass expansions 0-0, which are known 
to work very well in specific examples fH|,|l5]| . In other cases (e.g. the important 
two-loop QED corrections to muon decay |16|] and some electroweak corrections to 
B — ► X s j computed in |I7J), on ly gauge-invariant subsets of diagrams have been 
considered. Similarly to the case of the anomalous magnetic moment of the muon 
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considered in pi , the mixing B° — B° is a low-energy process which first occurs at 



one-loop level. Its two-loop electroweak corrections therefore require only one-loop 
renormalization and the relevant two-loop integrals can be expressed in terms of 
elementary functions, because the external momenta can be generally neglected. As 
will be illustrated in the following, this allows one to compute analytically the com- 
plete electroweak effects in a relatively simple way. The calculation presented here 
has therefore some interest in itself, independent of phenomenological applications. 
In particular, it exhibits the complete Higgs mass dependence and enables one to 
assess the validity of the heavy mass expansion in the case of electroweak boxes. 
Interestingly, our results show a very slow convergence of the heavy top mass ex- 
pansion in the case at hand. 



In addition to the strong motivation provided by the above three points, it should also 
be mentioned that our results represent an important subset of the electroweak corrections 
to the K° — K° mixing, whose treatment is more involved and will not be considered in 
the present work. 

This paper is organized as follows: in Sec. 2 we outline the strategy we have devised 
and the general features of the calculation; in Sec. 3 we describe the renormalization proce- 
dure in detail, giving explicit expressions for the renormalization constants we have used. 
Finally, Sec. 4 is devoted to a discussion of our results and of their main consequences, 
which are summarized in the Conclusions. 



2 Outline of the calculation 



2.1 The effective Hamiltonian 

The effective Hamiltonian TC^ =2 for the B° <-> B° transition can schematically be written 
as H 

where X t = V t * b V t d denotes the CKM factor, Q LL = by^(l — 7 5 )<i (g> hy^l — 75)0? is the 
relevant four-quark operator, and C B §({M}, fi) is the corresponding Wilson coefficient, 
to be calculated within perturbation theory. The latter contains all information on the 
heavy degrees of freedom, generically indicated by {M}, and depends additionally on the 
scale /i at which the operator Qll is renormalized. In principle, this scale can be chosen 
arbitrarily - physical quantities do not depend on it and therefore the /i-dependence of the 
Wilson coefficient will cancel against an analogous /i-dependence of the matrix element 
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(Qll) order by order in perturbation theory. It is physically natural, however, to choose 
\x of the order of since the b mass is clearly the relevant scale for the evaluation of the 
matrix element (Qll)- 

Including next-to-leading order (NLO) QCD corrections as well as two-loop elec- 
troweak corrections, the Wilson coefficient C BB assumes the form 

C BB ({M},fi b ) = r] 2B {S (w t ) + 5S ew } , (4) 

where the Inami-Lim function (w t = M^/Mf) 

1 - llw t + M 2 3\nw t . . 

S ° (ll "» = 4u,(u, - 1). + 2^31)3. (5) 



represents the leading order result [Tj|. The factor r\2B contains the complete NLO QCD 



corrections including the running of the Wilson coefficient down to the scale [i = m b H; 
f)2B is scheme dependent and its scheme dependence is canceled by analogous terms in 
the matrix element (see discussion in 0). In the NDR scheme and using an MS top mass 
evaluated at M t (see Sec. 3.1), f^B ~ 0.85. By 5S ew , instead, we denote the electroweak 
two- loop correction of order g 2 . The purpose of this paper is to calculate this quantity. 
For a better understanding of this calculation, and in order to introduce some relevant 
notation, we will first briefly recall some important features of the leading order (one-loop) 
calculation. 



2.2 Summary of one-loop results 

An important attribute of Wilson coefficients in general is their independence of the 
infrared region of the theory, i. e. "soft" physics. In particular, this means that one can 
calculate them using arbitrary kinematic configurations of the external particles (on-shell 
as well as off-shell), and arbitrary values for the light quark masses. In the case at hand, 
all quarks but the top can be considered light, and only the top has been integrated out 
in order to obtain the effective theory of Eq. (|3|). This freedom in the choice of the "soft" 
parameters crucially simplifies the calculation: in the following we will always adopt the 
simplest approach and work with zero external momenta and zero masses of the light 
quarks (except where these masses are needed as infrared regulator). 

Keeping this in mind, the one-loop amplitude for the process b + d — > b + d with 
vanishing external momenta and all light quark masses set to zero can be written as 

M lloop = t^t^t E AiA,- ^ Qll (6) 
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where Aj = V* b Vid and the sum is over the leading order box diagrams with internal quarks 
i and j = u,c,t), calculated in n dimensions. S^' is the general box function given 
e.g. in Eq. (2.2) of |J. Using m u = m c = and the unitarity of the CKM matrix, one 
obtains 



-i 9^1 
16vr 2 SMI 



M 



Hoop 



£(*>*) _ 2 S {t ' c) + S {c ' c) 



Qll- (7) 



The survival of the sole A t term in Eq. flj]) is a consequence of the hard (power-like) GIM 
cancellation mechanism which takes place in the case of electroweak boxes. For later 
convenience, we present here the explicit result for the function 

S(w t ,fi 2 /M 2 ) = <S (M) -2S^ + = S (w t ) + eS 1 (w t ,fl 2 /M 2 ) + 0(e 2 ) (8) 

where e = 2 — n/2 and we have included the 0(e) terms: So(w t ) has already been given 
in Eq. @, and the function Si reads 

_ 2 2 3 - 33w t -4 Wi 2 2 + 15^ t 31n 2 u;i 

Si{w t ,fl /M t ) = — 7T^- + T7 -rrlniOt- 



- l) 2 4(ty t - l) 3 4(w t - 1) 

M} 



+S (w t )ln-^, (9) 



with fl 2 = 47re" 7B /i 2 . 



Before proceeding further, some notation should be introduced. The electroweak cou- 
pling is generically denoted by g, the sine and cosine of the Weinberg angle by s and 
c, respectively. At the order we are working at, however, it is at some point necessary 
to specify the renormalization schemes (to be discussed in the following) in which these 
quantities are defined. Specifically, a hat (g, s, c) will always denote MS scheme quan- 
tities and the subscript W (g w , s w , c w ) the electroweak on-shell scheme. Furthermore, 
throughout the paper the following short-hand forms for the ratios of masses are used: 



Ml M 2 , M% 

w t = — % , zt = — I , h t = —| , (10) 
* M? ' * M t 2 * M? ' V ; 



where M w , M z , M t and M H are understood as on-shell masses. 



2.3 The matching procedure at two loops 

We are now ready to describe the strategy we have followed to perform our calculation. As 
a first step, we compute the complete electroweak corrections to the amplitude, setting 
all external masses and momenta and the internal light quark masses to zero. After 
renormalization, these corrections are of course ultraviolet finite, but some of the diagrams 
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containing the photon need the introduction of an infrared regulator. Following the QCD 
analysis of ||, we choose it to be a common internal light quark mass. The renormalized 
two-loop amplitude can then be expressed as 

- M2 '-=(l6^fl A,2,g - (U) 

where A^ is given by the sum of (i) the unrenormalized one-particle irreducible two- 
loop diagrams shown in Fig. 2a-h, (ii) the various counterterm contributions A^ (to be 
discussed in Sec. 3), (iii) a contribution Adt from the one-particle reducible diagrams 
depicted in Fig. 2i: 

A( 2 ) = AZ r oT + + a m w + + A^ + A c Tad + A DT . (12) 
The last contribution, Adt, can be gleaned from fl8| ; it originates from a finite but 



gauge-dependent subset of two- loop diagrams and was already considered in |J. We will 
discuss it in detail in Sec. 2.5. 

As a second step, we consider the effective theory of Eq. @, in which all degrees of 
freedom of the Standard Model (SM) with M > M w are decoupled, and perform the 
matching between the full theory and this effective theory at first order in QED. The 
amplitude in the effective theory is 



^TI^fl^M "' (13) 



whereQ 




A?//V) = \* ( ^ qt - In - 2 ) . (II ) 



Here, q t = m 2 q /M 2 (m q denotes the generic mass of the internal light quarks) plays the role 
of an infrared regulator. The electroweak correction to the relevant Wilson coefficient, 
proportional to the difference — A^ D (/i), is free from infrared divergences, as the 
dependence on q t cancels against analogous terms in A^ 2 \ The Wilson coefficient, however, 
depends crucially on the scale \x at which the operator Qll is renormalized. It is clear that 
our analysis also requires the knowledge of QED corrected matrix elements at the same 
scale [i. At least in principle, they could be provided by dedicated lattice calculations, 
but from a practical point of view it is likely that, since the QED corrections turn out to 
be very small, their impact would be irrelevant. 

2 In calculating Eq. (|lj) we have employed the projection method described in the next subsection, thus 
conforming to the same choice of evanescent operators as in || . A detailed discussion of the definition of 
evanescent operators can be found in j27l [28| . 
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We now consider the evolution of the Wilson coefficient down to the bottom mass 
scale. At lowest order in otQED, this amounts simply to setting fi = fib, with fi b pa m b 
in Eq. (|13|) . By adapting the available QCD calculations for the NLO anomalous 
dimension of the operator Qll to the case of QED, the evolution of the Wilson coefficient 
down to fib could in principle be performed consistently at NLO. However, as the QED 
coupling is very small and the QED logarithms have a negligible effect on the result, 
this simple approach (setting fi = fi b ) seems to us sufficient. Indeed, we notice that the 
second term in the parenthesis on the r.h.s. of Eq. (14), corresponding to the leading 



logarithmic QED correction, provides only a —0.3% correction to the one-loop result in 
the case fib = 4.8 GeV. 

Clearly, QCD corrections modify the evolution of the Wilson coefficient if the renor- 
malization group equation is solved keeping both QCD and QED effects into account. 
Using a <C a s , the formalism is well known [H], 20], and one finds that the second log- 



arithmic term in Eq. ([14]) should be multiplied by the QCD screening factor r}', given 
by 



1 1 \ (a s (M tJ 



TO 



V ' A)U(/i 6 ) a s (M t ))[a s (fib)J [ ln Mi) ] (15) 

the rest of the Wilson coefficient undergoes the usual NLO QCD evolution summarized 
by f]2B- Using (3q = y and the QCD anomalous dimension 7 = 4 ||, one obtains 
the numerical value rf ~ 0.9. This modest reduction of an already small effect can 
be approximately taken into account by noting that rj pa Tf?,B and including the QED 
logarithm In [xl/Mf {without the screening factor) inside the curly parenthesis of Eq. (|J). 
As a consequence of this discussion, we identify the electroweak correction SS ew of Eq. (f|) 
with the Wilson coefficient evaluated at a scale fib, which for definiteness we set equal to 



SS P . 



16tt 2 



A® - A«//V = m b )\ . (16) 



2.4 Calculation of the two-loop diagrams 



As far as the two-loop diagrams contributing to A^ in Eq. flTED are concerned, we need 
to consider all the topologies shown in Fig. 2, which correspond to a very large num- 
ber of diagrams. We have therefore decided to use the Mathematica |2l|] package 
FeynArts 1. 



.2 



p^] to generate automatically all the two-loop amplitudes. Dirac algebra 
operations, reduction to scalar integrals, and substitution of the scalar integrals have all 
been performed independently in two different ways involving various combinations of the 
Mathematica packages Tracer [23[] and ProcessDiagram |24]], and of Form ]p5 ]. The 
calculation has been performed in the 't Hooft-Feynman gauge, but we have checked the 
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£z independence. As we neglect all external momenta, the complete result for the un- 
renormalized two-loop amplitudes can be written in terms of two-loop vacuum integrals, 
which admit a relatively compact representation in terms of logarithms and dilogarithms 
of the internal masses (see e.g. [p6|). The final expression for the two-loop diagrams, 



however, is very long, due to the presence of four different heavy masses. Hence, we have 
decided to provide only approximate formulas that allow the reader to reproduce our nu- 
merical results with high accuracy. On the other hand, we feel that it may be more useful 
to explain in some detail the renormalization procedure that we have adopted. This is 
done in Sec. 3, where we stress the importance of keeping it as simple as possible in order 
to reduce the number of terms we have to deal with. All partial and final results are 
available in full analytic form from the authors. The complete analytic result has also 
allowed us to perform various heavy mass expansions, which will be presented in Sec. 4.. 

The computation of the contributions of individual two-loop diagrams to the Wilson 
coefficient of the operator Qll can be considerably simplified by the implementation of a 
suitable projection procedure. In our case, such a projection can be achieved by forming 
appropriate traces (see e.g. 2l| for a nice explanation of the method): writing a generic 



two-loop amplitude as 

M = bD 1 d®bD 2 d, (17) 

where D\ and D 2 are arbitrary Dirac structures with saturated Lorentz indices (for in- 
stance Di (g) D 2 = 7^7 A 7 P (1 — 75) ® 7^7a7 p (1 — 7s)), the projection on the operator Qll 
reads 

M - f 1 + tr [T M (! + 75)£>i7m(1 + 75)^2] Qll , (18) 



256 V 2 

where terms up to 0(e) have been included. Furthermore, since Qll represents the only 
operator contributing to the B° — B° mixing at the order we are working at, Eq. ( |T8"D is 
in fact an identity. 

A few additional comments are now in order: 



• We work in the framework of Naive Dimensional Regularization, a choice justified 
by the absence of any 75 ambiguity in the two-loop graphs we have considered. In 
particular, no closed fermion loop appears, apart from the ones in the self-energy 
insertions of the topologies in Fig. 2a, which pose no problem. 

• The number of diagrams is significantly reduced by the following simple observation: 
for a given topology, one can consider subsets of diagrams characterized by the 
number of light quark propagators they contain and by the lines to which the light 
quarks are assigned. Since all light quarks are treated as massless, diagrams within 
such subsets differ solely by the CKM factor, and can be grouped together, reducing 
the total number of diagrams to be actually computed. Using the unitarity of the 
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CKM matrix, the overall CKM factor turns out to be always proportional to Xf , as 
expected from the above-mentioned hard GIM cancellation. 



We perform the calculation in the electroweak scheme of Ref . [ 2§| , which uses MS 
couplings and on-shell masses as basic parameters. The transition to other popular 
schemes amounts to a finite renormalization of the electroweak coupling and is 
therefore straightforward: we discuss it in detail in the next section. The residual 
MS scale and scheme dependence allows us to gauge the remaining ambiguity due 
to the truncation of the perturbative series, as discussed in Sec. 4. 



2.5 Double triangle diagrams 

As mentioned above, the full electroweak two-loop correction receives also a contribution 
Ae>t from the one-particle reducible topologies depicted in Fig. 2i. This contribution is 
finite but gauge dependent. The two penguin diagrams may be connected by 7, Z, H, 
or <ft°, but only the bdZ vertex contributes to Aqt for vanishing external momenta and 
masses. Denoting this vertex by one finds in the 't Hooft-Feynman gauge 



<7 3 A 



^ = 7^^)^(1-75^, (19) 
with 

CW = 8^^-8K^ hm - <20) 
This result already includes the on-shell wave-function renormalization of the external 
legs. For the "double triangle" contribution Aqt we therefore obtain 

A DT = 16 C{w t f. (21) 

Numerically, g 2 /(lQir 2 ) x (A DT /S ) « +1.1%. 



3 Renormalization 



In renormalizing the two-loop amplitude our aim is to attain the maximal simplicity. 
We avoid all wave function renormalization of the internal lines, and choose a particu- 
larly simple procedure for the unphysical sector. All masses are defined on-shell and for 
the electroweak coupling we use the MS scheme 
connection to other schemes. 



29 1, although we explain in detail the 
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In the following we give explicit expressions for the various counterterms. They are 
written in terms of logarithms and of a single function B (x, y), which is defined through 
the one-loop integral 



-Re(g 2 e 7 ) e J 



d n k 



W 2 [k 2 - ml][(k - q) 2 - m|] 



B °(f>f)+°^ ( 22 ) 



whose analytic form is well known. The 0(e) part of the counterterms is not needed. To 
the reader's convenience we report here the explicit expression of Bq for the three special 
cases that are needed in our calculation: 



B (l,x) = 

B o (0,x) = 
B (1,0) = 

where the function a(x) is given by 



2 lnx a(x), 

2 2 K h 

2 — x In x — (1 — x) In 1 1 

2, 



a[x) 



2y/4x — x 2 arctan y 4/x — 1, < x < 4, 

x > 4. 



y/x 2 — 4 x In 



(23) 



3.1 Top Mass Counterterm 



The complete on-shell top mass counterterm is gauge-independent only after the inclusion 
of tadpoles, whose explicit expression can be found, for instance, in ||30|| . Here we report 
only the two-point function contribution to 5M t in the 't Hooft-Feynman gauge. It can 
be written as 



SM t 



9 



V ' 

M t 16vr 2 \hl? 



w t c z 



14 



0(e) 



where F t is given by 



F t = — 



4 — h t . . h t — 7 6 — h t , . , 
■a(/?, 4 ) H 1 — — ^ In h t 



16 



1 16 



1 



1 

w t 



16 

2^0(0,^) -^-±1(1- hi Wt ) 



1 

+ ^7 



+ 



Z t 



1 , /4 2 1\ 2 \ „ 1 s 2 4 



+ (2 + ^)(^-3j^5o(M t ) + - + y --, 4 



(24) 
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Replacing M t — > M t — 5M t in the LO result (pi), one finds the following top counterterm 

unren. 
2loop ' 



contribution to be added to A un " " 



with 5" = dS(w t , /i 2 /M^)/dw t , S' = dS /dw t , the function So(wt) as given in Eq. (0), 
and 5 defined by 

SK,^7<) = S{w t ,w t fx 2 /Ml), 

where S(w t ,p, 2 /M 2 ) has been given in Eq. (||). 

The top quark mass is therefore renormalized on-shell for what concerns electroweak 
effects?], while fj 2 B ~ 0.85 in Eq. (H) implies the use of a MS mass M t (M t ) as far as QCD 
effects are considered, in accordance to the standard convention ||. Using the LO QCD 
relation between on-shell and the MS mass and the pole mass value M t = 174 GeV JJJ, 
we find M t (M t ) = 166 GeV, which will be our input in the following and will be denoted 
for simplicity just by M t . 



3.2 Renormalization of the unphysical sector 

Before considering the renormalization of the W mass it is necessary to explain the pro- 
cedure we have followed for the unphysical scalars. As is well known, the renormalization 
of the unphysical sector is not independent from the way the physical sector is treated. 
Indeed, the renormalization procedure must respect the Slavnov- Taylor Identities (STI) 
which are induced by the local gauge invariance of the original Lagrangian before sponta- 
neous symmetry breaking. According to the organization of the calculation, it is possible 
to use different procedures that respect the STI's and are particularly convenient in or- 
der, for example, to minimize the number of counterterms to be considered. Of course, 
physical amplitudes are independent of the chosen procedure, and this can be used as an 
additional check of the calculation. For the problem at hand, the discussion can be kept 
at the one- loop level. 

If we split the unrenormalized W polarization tensor into transverse and longitudinal 
parts 

iCfe) = [aT - ^r) a ww (q 2 ) + ^n^( g 2 ), (26) 

3 A MS top mass renormalization in the electroweak sector would make the renormalized top mass 
dependent on the precise value of the Higgs mass, which is unknown; we do not consider it here. 
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and we denote by H w <t> the two-point function for the mixing between the W and its 
pseudo-Goldstone boson, 0, and by 11^ the self-energy of the latter, we obtain the fol- 
lowing STI in a general gauge (see for ex. |3T|.p2| ) 



q 2 (K, w (q 2 ) + 2M w U w<p (q 2 )) + M 2 W Tl^(q 2 ) + M 2 W T = , (27) 

where T represents the tadpole contribution. In particular, at q 2 = 0, the first two terms 
vanish, and the STI implies the cancelation between IL^,(0) and the tadpole contribution. 
This uncovers the connection between the renormalization of the Goldstone boson mass 
and the one of the tadpole. As anticipated, the counterterm contributions to the various 
terms in Eq. (|27|) must also respect the STI. In practice, the usual tadpole renormalization 
that minimizes the effective potential and consists in removing all tadpole graphs implies 
the subtraction of IT^O) from the two-point function of the pseudo-Goldstone bosonF]. As 
mentioned above, we adopt the physical mass to define the masses of the vector bosons. 
It is therefore convenient to renormalize the longitudinal component of the two-point 
function of the W in the same way as the transverse, using 5M^ = Re A WW (M^). For 
the other two-point functions different options are possible, which all respect the STI, 
and are equivalent at the level of physical amplitudes. They correspond to different ways 
of renormalizing the gauge-fixing parameters. 

One possibility consists in assigning no counterterm to the W — <fi transition; in the 
't Hooft-Feynman gauge this corresponds to renormalizing the masses of the vector boson 
and of the associated scalar boson in the same way {bare gauge fixing). Of course, the mass 
of the scalar boson will still need a supplementary subtraction at q 2 = 0, corresponding to 
the tadpole contribution. This choice clearly verifies Eq. (^) at q 2 = M^, leaving room 
for a further arbitrary wave-function renormalization, which we avoid altogether, as the 
W boson appears only inside the loops. Because of its simplicity, this is our preferred 
option: it amounts to fixing #M| = 8M^ + T and is the closest to the naive parameter 
renormalization. 

Another possibility, for example, would imply the subtraction of the first two terms 
of the Taylor expansion around q 2 = of the individual two-point functions in the 
external momentum, and it would obviously respect the identities, as the unrenormalized 
self-energies do. A counterterm for the W — <p transition is now involved. We have 
explicitly verified that all the renormalization options that satisfy the STI are equivalent 
at the level of physical amplitudes. 

Finally, we recall that the renormalization of the three point function du4> is fixed by 
the Ward identity that links the Yukawa coupling counterterm to the gauge coupling and 
fermion mass renormalization. 

4 This point is nicely explained in Taylor's book |53|, see section 14.6. 
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3.3 W mass counterterm 



Following the strategy outlined in the preceding subsection, we use the same counterterm 
to renormalize the W and the pseudo-Goldstone boson mass terms. In the 't Hooft- 
Feynman gauge, one has for the two-point contributions to the W mass counterterm (see 
e.g. 0): 



SMI 
Ml 



ReA ww {Ml) 



16tt 2 



M? 



1 



2w t 



(28) 



where the finite part Fw is given by 
(I 1 





1/1 h t \ 7 1 

^ 2+4 + — 



In fit 



-8c 2 



- 1 {-- 


h 2 A 
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wl I 2 
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6w t 
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h 2?- 


12c 4 



Replacing M w 
contribution 



M W —SM W in the LO result Eq. (0), one obtains the W mass counterterm 

(S-w t S') + (S -w t S^)F w \. (29) 




2w, 



1 + 7 

c 2 6 



Here again the notations S' = dS(w t , ft 2 / M^) / dw t and S' = dS /dw t are understood, and 
the functions So(wt), S(w t , p 2 /M 2 ) are given in Eqs. @ and (||). We treat the W mass as 
a fundamental input parameter, whose value is taken directly from the experiment: M w = 
80.385 ± 0.065 [0]. Although the determination of M w from M z , G^, and a is presently 
more precise than the experimental one after including all theoretical uncertainties |3^ 
this is certainly sufficient for our purposes. 



3.4 Tadpole contribution 

As mentioned above, the renormalization of the pseudo-Goldstone boson mass term can 
be performed (in the 't Hooft-Feynman gauge) in the same way as the one of the W mass, 
apart from an additional tadpole contribution, whose physical origin has been already 
explained. In practice, this means that the one-loop diagrams containing scalars induce 
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a further counterterm contribution that can be easily calculated from the complete one- 
loop tadpole (see for instance Eq. (12) of |30[| ) and the one-loop amplitudes involving only 
scalars. Here we just give the final result for this additional contribution, again in the 
't Hooft-Feynman gauge: 



A c 




E, 



Tad 



e 



Tad 1 



Ol 



(30) 



with 



Ti 



Tad 



Tad 



-1 + 19w t 1 - 6w t 



4K - 1) 3 

-13 + 43w t 



2(w t - l) 4 
4 - 19w? 



Aw 2 

lnw t , 



6wt — 4w 2 




In w t + T \n^ 



Mr 



]c 2 w t 



\nz t . 



3.5 Coupling counterterm 



The SU(2) coupling counterterm in the MS scheme is (up to ln47r and je constants) 

Sg _ g 2 19 

g 16tt 2 12e ' 

Performing the renormalization go = g — 5g, we obtain the following contribution to A^: 

19 

Al = --S{w t ,^/Mf). (31) 

In the MS scheme g(M z ) can be calculated from a{M z ) = (128. 1)" 1 and s 2 {M z ) w 
sm 2 9 l e J} - 1 x 10" 4 = 0.23145 |3f, obtaining g 2 (M z ) = 0.423842; this is the value that 



will be used in the following numerical calculations. 

Let us now consider what happens in different renormalization schemes. In any scheme 
X but MS the counterterm would have a finite part , namely 

5g x _ g 2 19 F x 

g 16vr 2 12e 9 ' 

so that the difference between the result in a scheme X and the result in the MS scheme 
at the scale fi g = M z is just 

6S x = -4S (w t )Ff, (32) 
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which should be added to 5S^, also evaluated at \i g = M z . It is well-known that the 
use of the electromagnetic fine structure constant a(0) = 1/137.036 to normalize the 



electroweak coupling would introduce large and uncompensated mass singularities [36] in 
the two-loop result. Therefore, the coupling g must be normalized at short distances. 
This is the case, for instance, if we use G M and M w : g 2 = 4 a/2 M^. Then we find 



F G. A + 



9 



(33) 

MS/ 

where Af w is a function of the top, Higgs, W and Z masses, as well as of the couplings, 
which summarizes the electroweak corrections to the muon decay in the MS scheme and 



is given explicitly in [29]. The second term in parenthesis, instead, represents the finite 
part of the conventional electric charge counterterm fl36l , evaluated at the MS scale fj, g . 

depends mildly (logarithmically) on M H and M t \ numerically it is very small, as can 
be seen comparing g 2 obtained from G^M^ and the MS value of g 2 (M z ): they differ by 
half a percent. 

We will consider later also another possibility to define g, namely through the relation 
g 2 = 47r a(M z )/ s^. Here we have adopted the notation = 1 — M^/M 2 for the on 
shell definition of the weak mixing angle |36[ and a(M z ) stands for the MS running 



electromagnetic coupling at the scale M z .| In this case 

1 c 2 

2 s 



^oTfAp, (34) 



where A/5 is another function of the heavy masses given in Ref. ||29|| . Unlike Eq. (|33|), 
however, Fg° s is very sensitive to M t , due to its quadratic dependence enhanced by the 
factor c^/ s^y. Indeed, the MS and on-shell definitions of sin 2 9 W , and so the corresponding 
g 2 , differ by almost 4% ! We can therefore anticipate that this normalization choice will 
introduce large two-loop corrections when compared to the MS scheme. Moreover, as the 
origin of the large counterterm contributions has nothing to do with the process at hand, 
we can also expect this choice to give unnaturally large two-loop corrections. 



3.6 Wave function counterterm 



The left-handed quark fields of the down type (i. e. the external fields) are renormalized 
according to 

df, (35) 

5 To keep the notation simple, we use here the MS running electromagnetic coupling a(M z ) ~ (128. 1) -1 
instead of the conventional a(M z ) ~ (128.9) -1 ||, which is defined by subtraction of the renormalized 
light fermion photon correlators only, and is more generally used in the context of the on-shell scheme. 
The difference is not negligible, and must be taken into account in the two-loop corrections. 



(dl 



(Z L d 



,V2 
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and the (matrix-valued) wave-function renormalization constant = 1 + 5Z% is de- 
termined as described by Aoki et al. [31|. Denoting the sum of all one-loop diagrams 
contributing to the transition di 
(L,R=(l T l,)/2) 



dj with external momentum p by iEy(p), and writing 



EM = EUp 2 )i>L + EfAp^R + Sj/IKI + mjR) 



(36) 



this prescription implies in the limiting case of vanishing masses of the down-type quarks 
the customary relation 

{5Z L d h = -S§(0). (37) 

In this equation we have considered only the hermitian part of the wave function renor- 
malization constant. For what concerns the antihermitian part of 5Z%, it is not defined 
in the limit of vanishing quark masses that we have adopted from the beginning. One 
can use the renormalization prescription of Denner and Sack |37] for the CKM matrix 
elements and remove it completely. In the case of the up quarks, the top mass cannot 
be set to zero, but the antihermitian part of the wave function renormalization vanishes 
because of the GIM mechanism. 



By calculating the relevant diagrams, we obtain in this way 



where 



9 / H 



16tt 2 



_ 2 \ £ 



(3* 



A = — 



B 



e 136 18 J 72 12 18 V36 9 9 1 



+ gS 2 \nq t , 
1 1 



3(1 + w t ) w t -A 

+ — l n Wf 



eAw t 8w t (l-w t ) 4(1 - w t f 



(39) 
(40) 



The corresponding counterterm contribution AS, is obtained as follows: the amputated 

amplitude M.xi oop for the transition k + d—>b + d(k = d,s,b) is in the notation of Sec. 
" ™ given by 



M (k) -i 9 
Hoop 167r 2 8M 2 



(41) 



w i,j 



with X\ = V* k Vid- Using the unitarity of the CKM matrix, this simplifies to 



L 

-i g 



M (k) 

Hoop 167r 2g M 2 



[Xfht [S {t ' f) - 2S {t ' c) + 5 (c ' c) ] + 5 kd X t [S^ c) - 5 (c ' c) ] } . (42) 
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(k) 

The second term contributes only for k = d, and it originates in the unitarity relation X t + 
X'f' + \w = Sf-d- The counterterm contribution due to the wave function renormalization 
of the external 6-field is then 



bk JVl Uoop- 



(43) 



Combining (RSI) and one thus obtains 



ci 1 / 



^ 2 

with S(w t ,p> 2 /M 2 ) as given in (||), and 
U(w t ,p, 2 /M 2 ) = S^-S^ 
+ e 



(A + B) S(w t , fi 2 /M 2 ) + B U(w t , tf/M 2 ), 



1 



+ 



1 w t In w t 

w t -l (w t - l) 2 
2-w t 



1 + e In 



2(w t - 1) 2(w t - l) 2 



lnw t + 



w t In w t 



M? 



+ 0(e 5 



(44) 



An analogous reasoning holds for the wave function renormalization of the other three 
external fields, and the final counterterm contribution is simply 



4A 



c,l 



(45) 



4 Results and discussion 



Combining the counterterm contributions of Eqs. (p5|) , (^-R_Tj), (|^) with Eq. (|2l|) and 
with the unrenormalized two-loop contributions yields our final result in the MS scheme. 
The size of the effect can be seen in Fig. 3, where the following input values have been 
used 0: M w = 80.385 GeV, \i g = M z = 91.187 GeV, M t = 166 GeV, (i b = 4.8 GeV, 
s 2 (M z ) = 0.23145. These are the values we will use throughout this section. The elec- 
troweak correction is less than 1% of the LO result for any value of the Higgs mass below 
1 TeV, and is particularly small for a light Higgs, namely in the region preferred by recent 
global fits |f7l,[35|. The logarithmic asymptotic dependence for large Higgs masses is al- 
ready evident around 1 TeV. In Sec. 4.3 we will present a simple approximate formula that 
reproduces the complete result very accurately. Before, however, we consider the heavy 
mass expansions of the full result and study the residual scheme and scale dependence. 
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Figure 3: The electroweak correction 5S ew in the MS scheme, normalized to the L0 result, 
as a function of the Higgs mass. Inputs as in the text. 



4.1 Heavy mass expansions 



We first consider the heavy Higgs expansion. In the case of a very heavy Higgs boson, 
M H 3> M w ,M t , we have verified that the leading behaviour is logarithmic in M H . In 
particular, the MS scheme result reads 



6S: 



hh 



fr 



r l n 2 ht 

32 



167T 2 W 2 



+ 



27w? + 30w t 



48K - 1)' 



5 w\ (3w t — 5) lnw t 



16 (w t - 1) 



In ht + 



(46) 



where the ellipses stand for constants or terms suppressed by inverse powers of the Higgs 
mass. Of course, in another scheme X the Higgs dependence contained in 5Sx will 
change the In ht term in the above expression. The absence of a leading quadratic Higgs 
behaviour in the complete result can be explained in the context of an effective theory 
approach where the Higgs boson is integrated out . 



If we are interested in the leading heavy-top corrections, 0(g 2 Mf/M^), we may work 
in the framework of a Yukawa Lagrangian where the heavy fermions couple only to the 
Higgs boson and to the longitudinal components of the gauge bosons, a situation which 
corresponds to the gaugeless limit of the SM 111] . In the covariant gauges, it suffices to 



consider only the diagrams that involve exchanges of scalar bosons. Even in the case of 
a heavy Higgs boson, in which the scalar coupling A ~ g 2 M 2 /M^ cannot be neglected 
with respect to the Yukawa coupling of the top, only a few tens of two-loop diagrams 
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contribute. Notice that in this framework the mass of the Higgs boson is completely 
arbitrary. 



Computing the various counterterms from the scalar one-loop diagrams only (which 
corresponds to the heavy top limit of our renormalization constants), adding these to the 
sum of the scalar two-loop diagrams, and expanding in small Wt, we get 



SS. 



lead 



r 



167T 2 Wf 



-44 + 112h t - 66/i 2 + Uh 3 t - h\ , f h\ 52 - 16 h t -9h 2 t 
64 \J) + 32 

22h 2 - 20/i? + hf 



(ft, - 1)< (4 - 14 VhftB ht) _ — 13 h , hi 



16 

h t (U + 9h t ) 



32 



8 + 22h - 20^ + hj 2 ,„ 
In ht 7T + 



64 
/k-4 



192 



32 



a(h t 



(47) 



Here we have indicated the dilogarithmic function as Li 2 (x) = — Jq dt ln ^ t - , 
duced 



4JiZi Cl 2 (2arcsin x /i), 



-4Li 9 



+ 2 In 



-In 2 (4,2) + tt 2 /3 



0<z<l, 

z > 1, 



and intro- 



(48) 



where A = y 1 — \ and Cl 2 (x) = ImLi 2 (e ia: ) is the Clausen function. It is a welcome 
check of our result that Eq. (|47|) coincides with the limit w t — > of the complete result, 
which is scheme and scale independent. The limit for M H — > of the square parenthesis 
in Eq. (f47|) is y — One can also verify that the heavy Higgs limit of Eq. ( f4"7] ) coincides 
with the heavy top limit of Eq. (f46|) . 



Numerically, 5S l ™ d /S = 1.12% (for M H = 100 GeV), which should be compared to 
0.18% obtained using the full result. Clearly, the leading term of a heavy top expansion 
does not approximate well the complete result. This is not surprising since, for the 
measured value of the top quark mass, the heavy top expansion of the LO result, 



So(w t ) 



1 



4:W t 



9 
4 



f 15 9 , \ ^/ 
mw t - w t I — + - In w t J + U{w, 



(49) 



must be pushed up to the fourth term in order to reach a 10% accuracy. An example of 
slow convergence of the heavy top expansion can also be found in the 0(aa s ) corrections 
to the Z —>■ bb decay amplitude |39|], although in that case the leading top contribution 
dominates, because of large cancellations among the subleading terms. 

It is therefore interesting to investigate further the convergence of the heavy top ex- 
pansion in a case where the complete two-loop contributions are available in analytical 
form. We have expanded our complete MS result in inverse powers of the heavy top mass, 
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Figure 4: (a) Renormalized box contributions to SS ew in the MS scheme, normalized to the 
LO result, as a function of the Higgs mass. The solid line represents the complete result, 
the dotted line the leading top approximation, and the dashed (dashed-dotted) line the heavy 
(light) Higgs top expansions up to 0(1/M t ). Inputs as in the text, (b) The same, but for 
M t = 350 GeV, demonstrating the much better convergence of the heavy top expansions in 
this case. 



up to the third term, namely excluding only contributions which are formally 0(M w /M t ). 
With respect to the Higgs boson mass, we either consider it much smaller than the top 
mass, or comparable, i.e. heavy compared to M w . Consequently, we obtain two expan- 
sions, one valid for small M H , the other for large M H . In Fig. 4, we compare them with 
the leading top result of Eq. (|4~7D and with the complete result; we display the various 
results as functions of M H for the two cases M t = 166 GeV and M t = 350 GeV, respec- 
tively. We show only the contributions of the renormalized box diagrams, i.e. we leave 
aside the double triangle contribution Ajjt', the latter does not depend on M H and its 
top expansion converges also very slowly, as can be seen from its explicit expression of 
Eq. (|2T1). We observe that in both cases the light and heavy Higgs expansions match quite 



well around M H « 60 and 100 GeV (similar to what happens in [12]), but start to converge 
to the complete result only for a very heavy top, with M t = 350 GeV. For even heavier 
top masses, the expansion up to 0(1/M t ) terms provides an excellent approximation. 

The slow convergence of the heavy top expansion of the two-loop contributions may 
be related to the slow convergence of the expansion of the LO result. It is therefore quite 
different from what happens in the case of the calculation of precision observables [ T% TA 



where the top contributions originate from two-point functions. The heavy top mass 
expansion for the latter seems to work remarkably well up to two-loop |15j and has been 



checked even at the three-loop level in the case of mixed 0(aa s ) corrections . 
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Figure 5: Dependence on the MS scale of the LO and of the sum of LO and two-loop results 
for M H = lOOGeV. Inputs as in the text. 

4.2 Scale and scheme dependence 

In the MS scheme, we have so far fixed fj, g = M z . However, by varying \i g around M w , we 
can compare the dependence on the renormalization scale n g of the LO result with the one 
of the sum of one and two-loop results. We recall that the LO result is expected to exhibit a 
strong scale dependence because of the factor g 4 (/z s ). Our two-loop correction cancels the 
leading logarithmic dependence on /i g , and the remaining much weaker scale dependence 
can be interpreted as an indication of the importance of higher order effects. This is 
illustrated in Fig. 5, where we show that in varying the scale fi g between M w /2 and 2M W 
the combination g 4 So(w t ) varies by about 5%, while the combination g 4 (So(wt) + SS ew ) 
undergoes a maximum variation of less than 0.1%. This dramatic reduction of the scale 
dependence demonstrates the importance of a complete two-loop electroweak analysis in 
order to reduce the normalization ambiguities of the LO result. 

Following the discussion of Sec. 3 on the normalization of the electroweak coupling g, 
we can also compare the LO result g A S (w t ) and the electroweak corrected g A (S (w t ) + 
5S ew ) in different renormalization schemes. Table 1 summarizes the scheme dependence 
before and after the inclusion of our two-loop calculation by comparing the three repre- 
sentative cases introduced in Sec. 3: MS scheme at the scale M z , g expressed in terms of 
G^M^, g expressed in terms of the on-shell sine and of a(M z ). The input values that 
we use are to a good approximation compatible with the precise calculation summarized by 
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scheme 


g 4 S 


g 4 (S + 6 Sen) 


MS 
OS 


0.4270 
0.4320 
0.4603 


0.4277 
0.4276 
0.4270 



Table 1: Scheme dependence: comparison between LO and electroweak corrected results for 
M„ = lOOGeV. 



the numerical formulas provided in Ref. []35| in the case M H = 100 GeV and M t = 174 GeV 
(equivalent to M t = 166 GeV). For our purposes it is therefore sufficient to evaluate nu- 
merically the extra terms 5S<js and SSq inverting the definitions of A/3 and Ar w : one 
finds F° s = (1 - s 2 w /s\M z ))/2 and if * = 1/2 (1 - 7ta(M z )/(V2 s 2 (M z ) G M )), from 
which 5Sos an d $Sgu can be computed using also G M = 1.16639 x 10~ 5 GeV~ 2 . Here 
we have implicitly assumed that the theoretical determination of the relations among the 
electroweak couplings is ideally accurate; although this is obviously not the case, the the- 
oretical errors involved are of the order of a few parts in 10 4 [[HJ and can be neglected at 
this stage. 

As could be expected after the remarks made in Sec. 3 about the scheme dependence, 
the difference between MS scheme and normalization is only 1.1% at LO, while the 
one between the MS and on-shell schemes is unnaturally very large for a purely elec- 
troweak correction at low-energies, almost 8%. After the inclusion of our two-loop correc- 
tions, the situation changes drastically. The maximum differences are now about 0.2%, 
which agrees with a rough estimate of the residual uncertainty, obtained by squaring the 
relative difference between s% and s 2 (M z ), that g 0.002. 

Because the on-shell scheme induces large and theoretically well-controlled radiative 
corrections 5S ew , however, it cannot be used to estimate higher order effects when the 
couplings are normalized in a different way. We have therefore extended the scheme 
comparison of Table 1 to different values of M H , limiting ourselves to the case where the 
result is expressed in terms of MS couplings or of G M . To this end, one needs to recalculate 
s 2 (M z ) for each M H from G M and M w , and to use the complete expression for F^". We 
find excellent agreement between the two schemes: the residual ambiguity is always well 
below 0.1%. 
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4.3 Approximate formulae 



In this subsection we present some very simple approximate formulae that reproduce with 
excellent accuracy the result of our calculation when the LO result is normalized to the 
MS coupling g(ng) and to G^M^. Unlike the original very long analytic formulae, they 
are suitable for a simple implementation in numerical analyses. 



In the MS scheme, we can rewrite 8S ew isolating the dependence of the complete result 
on the two scales fi g , at which the coupling g is normalized, and /x& ~ m&, at which the 
Wilson coefficient of the underlying effective theory is evaluated (C(w t ) is the penguin 
function given in Eq. (p0|)): 



SS. 



MS 



g- 



16tt 2 



19 , ft 



S (w t ) [ In r*. + 1§ 2 In ^ ) + A(M H , M t ) + 16 C(w t ) 



M 



M 



(50) 



The following simple approximate formula reproduces with high accuracy the analytic 
result for A(M H , M t ), which cannot be reduced to a compact form: 



A(M H ,M t 



12.20 + 52.74 



1 



0.977 In 



Mh 

100 



1.709 In' 



100' 



(51) 



where we have expressed the top and Higgs masses in GeV. For 80 < M H < 650 GeV, M t 
within la from the present experimental value, M t = 166 ± 5 GeV, and using s 2 (M z ) = 
0.23145 and M w = 80.385 GeV, Eq. (|51|) does not deviate from the analytic result by more 
than 1% which induces in the complete result a very small O(10 -4 ) relative error. 



In the case the LO result is expressed in terms of G^M^, only the fib scale dependence 
is left and the Higgs and top dependence is different due to the additional term SSq „, as 
shown in Fig. 6. The result can then be written as 



6S g, = V2G,M 2 W 



4tt 2 



\ S (w t ) In J* + B(M H ,M t ) + 16 C{w t f 



where 



B(M H ,M t 



-15.49-31.65 



166 



100 



100 



(52) 



2.296 In + 1.868 In 2 tt^, (53) 



and again M t and M H are expressed in GeV. In the same ranges of M t and M H that 
we considered for the function A(M H , M t ), B(M H , M t ) approximates the complete result 
with analogous accuracy. 
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Figure 6: The electroweak correction 5S ew when the L0 result is expressed in terms of G^, 
normalized to the LO result, as a function of the Higgs mass. Inputs as in Fig. 3. 

5 Conclusions 

In this paper we have computed the complete electroweak effects in the B° — B° mixing 
in the SM. We have used an effective theory approach in treating the QED corrections 
and neglected QED effects in the matrix elements. The calculation has been performed 
analytically, and the results expressed for arbitrary Higgs mass in terms of the very simple 
approximate formulae of Eqs. (|50|- |53|) . The discussion of Sec. 4 can be summarized in the 
following way: 

1. Unless the renormalization procedure introduces unnaturally large correction un- 
related to the process, the two- loop electroweak corrections are small, 0(1%), for 
any realistic value of the Higgs mass. This is of the same order of magnitude of the 
NLO QCD corrections. If one uses an on-shell definition of the electroweak coupling, 
however, they reach almost 8% for a light Higgs boson. 

2. The electroweak scheme and scale dependence of the result, which we have shown 
to be quite large at LO, is consistently reduced to the permille level. If G^M^ or 
MS couplings are used to normalize the result, the ambiguity is below 0.1%, much 
less than the present perturbative QCD error. 

3. The heavy top expansion converges very slowly in the case of the B° — B° mixing. 
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This was known at LO and we have verified the same pattern at the two-loop level. 
The fact that the leading heavy top limit does not approximate well the complete 
result suggests that it cannot be used as a reliable estimate of the full electroweak ef- 
fects in cases where already at LO it does not work properly. In particular, this 
refers to some weak decays for which the leading heavy top corrections have recently 
been computed: see [[40| (for the decay B — > X s ^) and || (for a class of rare decays 
including K — > -nvv). 



If one expresses the LO result in terms of and M w , as it is customarily done 0, the 
correction as a function of the Higgs mass is shown in Fig. 6 and the effect on the extraction 
of \Vtd\ from AM B<J is tiny. Indeed, \V t d\ ~ (So + 8S ew )~ 1 ^ 2 , so the electroweak corrections 
induce approximately a +0.5% change in the extracted value of \V t d\, for values of the 
Higgs mass below 400 GeV. 
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